PARTIAL-FRACTION DECOMPOSITIONS 
AND HARMONIC NUMBER IDENTITIES 



CHU WENCHANG 



Abstract. By means of partial fraction method, we investigate the decom- 
position of rational functions. Several striking identities on harmonic numbers 
and generalized Apery numbers will be established, including the binomial- 
harmonic number identity associated with Beukers' conjecture on Apery num- 
bers. 



1. Introduction 



The generalized harmonic numbers are defined to be partial sums of the Riemann- 
Zeta series: 



i/<">=0 and i/(™>=^i- for m, n=l,2,.... (1) 

k=l 

When m = 1, they reduce to the classical ones, shortened as _ff„ = Hyi\ 



If the shifted factorial is defined by 

(c)o = 1 and (c)„ = c(c 1) • • • (c + n - 1) for ri = 1,2,- •• (2) 

then we can establish, by means of the standard partial- fraction decompositions, 
the following algebraic identities: 



{x) 



ri+l 



fc=0 



3 



Multiplying H5al5b|l across by x and then letting a; ^ oo, we recover one identity 
among the hardest challenges claimed in W, Eq 16], [S! Eq 12] and |1(J[ Eq20]: 

n o 

Y.^-lf{l) {3(if,-i7„_,)^ + «+i7f2,)} = 0. (6) 

fc=0 



Date: February 1, 2008, The work carried out during my visit to Center for Combinatorics 
(LPMC), Nankai University (2005). 

2000 Mathematics Subject Classification. Primary 05A10, Secondary 11B75. 

Key words and phrases. Apery number, Binomial coefficient. Generalized harmonic number. 
Shifted rising factorial. Partial fraction decomposition. 



2 



CHU Wenchang: Partial-Fraction Decompositions and Harmonic Number Identities 



This has best exemplified tlie power of partial fraction method. For more general 
rational functions, we will invc;stigate their partial fraction decompositions in the 
second section, which involve the complete Bell polynomials (or cyclic indicators 
of symmetric groups) on the generalized harmonic numbers. Several further ex- 
amples and miscellaneous formulae will be collected in the third and last section. 
In order to facilitate consultation for readers, three short tables of the complete 
Bell polynomials on the generalized harmonic numbers will be presented in the 
appendices. 



2. Partial Fraction Decompositions 

For two natural numbers n and k with < fc < n, define two functions related to 
harmonic numbers by 

n ^ 

^^(^)^=E77— =^ Hii-k)=Hi%,-Hi'^ (7a) 

L=l ^'^ ^> 

n ^ 

W^W^=E(7T^ ^ w,(-fc)=i?i2fe + (-i)'<^ (7b) 

They come respectively from the logarithmic derivatives of the binomial coefficients 
hi \ (1 ~ (n-x\ , 



n\ X {x + k) 



Let a{tj be the set of partitions of I represented by ^-tuples of nonnegative integers 
(mi,m2, • • • ,m^) such that X^^^i kruk = i- Its subset of ^-partitions into m parts 
with Yfk=i "^fe = mis denoted by crm{i)- 

Theorem (Partial fraction decomposition). Let A, /i and n he three natural num- 
bers with X + {X — > 0. Then there holds the algebraic identity: 

(n!)^-^! - .x)(^ _^ kx{nY(n + kY^ ne{X,fi,-k) 
to ^ ^ ^ £^^!(x + fc)^-^ 

with the Q,-coefficients being determined by the Bell polynomials (or the cyclic in- 
dicators of symmetric groups) ; 

^ ^i^{xn.{x)-{-iynHi{x)Y' 
ne{x,^.,x) = (-i)^^!En-^ — (10) 

cr{e)i=l 

where the multiple sum runs over a{t), the set of £ -partitions represented by £ -tuples 
of nonnegative integers (mi, m2, • • • , mi) such that X]i=i = 



In particular, the 0-coefScients read explicitly as 

ne{X,n,-k) = £lJ2l[^ (11) 

a{e)i=l 
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Proof. By means of partial fraction decomposition, we can formally write 

where the coefficients C{k,£) are to be determined. Letting Vx = stand for the 
derivative operator with respect to x and then noting that 



we first demonstrate that for < £ < A there holds: 

Cik,i) = h\-k)h^{-k) X M^iifl^ (12) 
where the J7-coefficients are given by the following logarithmic derivatives: 



,^h^{x)h^'{x)^ 



^ h\x)h^{x) ■ 



For £ = 0, we have obviously rJo(A, /x, x) = 1 and that 



C(fc,0) = lim ft^(a;)/?,^(a;) X f7o(A,^,a;) 

x^ — k 



= h''i-k)h''{-k) X f7o(A, -k). 
Next for £ = 1, we can check (|12I13|I through L'Hospital's rule that 

- -.'-^>-'{^-|^: 

h^{x)hf'{x) -C(k.O) 
= lim 

= \im^Vx\^h^{x)ht'ix)^ 

= h^{-k)h>'{-k) X f7i(A, At, -k). 



Supposing now the truth of H12I13|I for ^ = 0, 1, • • • , m — 1 with m < A, then we 
have to verify it also for £ — m. Applying again the L'Hospital rule for m-times, 
we can determine the coefficient 



- m — 1 

^^[h\x)h>^{x)] 



- lim h\x)h^'{x) . . 

x^-k to! f^{x) 

= h\-k)h^^{-k) X ^rn{\fi,-k) ^ 

to! 

Based on the induction principle, we have confirmed that the coefficients in partial 
fraction decomposition are determined by p2ll3|l . 
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To complete the proof of the theorem, it remams to show that these coefficients can 
be calculated exphcitly through equation Hl()() and therefore ()ll(l (Bell polynomials 
and/or the cyclic indicators of symmetric groups). 



Manipulating the differential operation 

Vl+'^l^h^{x)h''{x)} V^J^h^{x)h''{x)} Vi^^h^{x)h^'{x)} 
h^{x)hf^{x) ^ h^{x)hf^{x) h^{x)hf^{x) 
we can derive for (|13|1 the recurrence relation 

ni+e{X,n,x) = ^^V,-XHiix)- ^lHl{x)}ne{X,^i,x). (14) 

It is trivial to see that fli{X,fj,,x) defined by (|10|l admits the initial condition 
ilQ{X,ii,x) = 1. If we can check that il.i{X,^,x) defined by (|10|l satisfies the same 
recurrence relation p4l) . then the validity of (|10|l would be confirmed for all the 
natural numbers i. 



Now substituting the RHS of lfTn|l into the RHS of (fTljl and then noticing the 
differential relations 

V^Hj{x) = +jHj+i{x) 
V^TLjix) = -jHj+i{x) 

we get the following expression 

i- ly+'elf\xn^ix) + nH,ix)]J2Y[ 



{xn^{x)-{-ly^iH,{x)j 



(. i=l 

<^ Xn^(x) - (-lYnHJx) } i . ■ 



mji™' X!-?"^. 

a{l) i=l i=l 



1 \Hj(x)-(-iy tiH,(x) 



(15a) 
(15b) 



In accordance with the combinatorial structure, each ^-partition enumerated by 
a{t) becomes a (1 + ^)-partition with a "j"-part being shifted to a "1 + j"-part for 
< j < ^- Vice versa, every (1 + £)-partition enumerated by cr(l + 1) reduces to a 
^-partition with a "1 +/'-part being replaced by a "/'-part for < j < ^. Then 
the sum over partitions should be reformulated accordingly. 



First, the line H15a|) with an extra part "1" yields a new factor Mi :— mi + 1. 
Then if — 0, for each j corresponding to the shift from part "j" to part "1 + j" 
displayed in line Ijl5b|l . the coefficient jmj is replaced by (1 + j)ilfj+i under two 
index substitution Mj := rrij — 1 and Mi^j :— mi+j + 1 for I < j < £■ Lastly if 
mi = 1, the coefficient £mi will be replaced by (1 + i)Mi^i with two summation 
index being substituted by Mi := me — 1 and Afi+f := mg. Summing up, we may 
combine H15a|l with l)15b(l and obtain the following expression 

According to 1)10(1 . the last expression becomes fli+i{X, ii,x), i.e., the left member 
of (O thanks for the (1 + £)-partition 1 + £ = Y.]tljMj. This confirms that 
the RHS of (|10|) satisfies indeed the recurrence relation H14|) . We therefore have 
established equality (|10|l and Hll|) . This completes the proof of the theorem. □ 
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In the theorem, multiplying the partial fraction decomposition by x and then letting 
x — > c», we derive the following harmonic number identity. 

Corollary (Harmonic number identity). Let A, /i and n be three natural numbers 
with A + (A — /i)ri > 1. Then there holds the algebraic identity: 

Ei-iriiyct^y^x-ii^i^.-k) = (16) 

where the ri-coefficients are given by the Bell polynomials (|10(l . 
If defining further two sequences by 

= n,ixM = = eij: (-i)^A™n^ (17a) 

^ ' a(t) i=l 

and then applying the Leibniz rule to H13() . we find the following convolution for- 
mula: 



£ 

il.e{X,n,x) = ^^(^ ^^■!u,{\,x)uji-i,{p,x). (18) 

t.=0 



Putting a; — —k, we write down the corresponding relation as follows: 

I 



■{X,fi,~k) - ^(;)n7,(A,-fc)u;,_,(M,-fc) (19) 



where w and uj are explicitly provided by the following formulae: 

w,{\-k) = A,0,-fc = £!VA™r[i — ^ 20a 



t{1) i=l 



u;e{fi,-k) = n,iO,fi,-k) ^ ■ (20b) 



3. Examples: Harmonic Number Identities 

By means of the theorem and the corollary, we will display several examples of 
partial fraction decompositions and the corresponding harmonic number identities. 

Example 1 (A = 1). For (5 = 0, 1, there hold partial fraction expansions 
I ,M-,5 (1 -^)^ _ ^( n\(n + k\^{-lf 

and the two corresponding harmonic number identities: 

n 



k=0 



i^fn\fn + k 
k 



6 



CHU Wcnchang: Partial-Fraction Decompositions and Harmonic Number Identities 



When fi — 0, the corresponding partial fraction expansion reads as the formula 
displayed in 0. For an alternative derivation, refer to the recent paper When 
11=1, the last binomial identity is a special case of the Chu-Vandermonde convo- 
lution formula on binomial coefficients: 



[ n ) ~ fc J V n- 



fc=0 



Example 2 (A = 2). For /i = 0, 1, 2, there hold partial fraction expansions 



E 



fe=0 



{x + fc)2 

X + k I 

and the three corresponding harmonic number identities: 



fe=0 



k 



For /i = 0, the corresponding partial fraction expansion has been given by 
Example 3 (A = 3). For /i = 0, 1, 2, 3, there hold partial fraction expansions 

"^^^l k) [ k ) \(x + fc)3 



+ 



(a; + fc)2 
1/2 



(3 + ^i)Hk - 'iHn-k - l^H„ 



{(3 + M)<+3ffr-^,-^i/r^J 



(2) 



r(2) 



anrf the four corresponding harmonic number identities 

n „ 



3/n + fc\A'j [(3 + n)Hk - 3Hn-k - l-l-Hn+k] ' 



k=0 



2) 



r(2) 



When /i = 0, the corresponding partial fraction decomposition and harmonic num- 
ber identity have been exhibited respectively in H5al5b|) and 10 • 

Example 4 (A = 4). For fi ^ 0,1, 2, 3, 4, there hold partial fraction expansions 



(^)n+l 



E 

fe=0 

1/2 

(x-|-fc)2 

1/6 

x+k 



k ) I {x + fc)4 



+3{(4+Ai)Hfc-4ff„_fc-MH„+fc}x{(4+M)-f^f *+4Jfi'ifc-Mffl^H^fc} 



and i/ie /jwe corresponding harmonic number identities: 

Q _ -A / n \4 /n + f [(4+M)Hfc-4H„_fc-pH„+fc]%2[(4+M)Hf >-4Hf/,.-pfff_^J 

^Vfcy V A: y \+3[(4+/^)_Hfc-4_H„_fc-/i_H„+fc]x[(4+M)fl"f*+4Hi!!fc-Mffi+fc] 
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Example 5 (A 



(n!) 



5). For /J, = 0, 1, 2, 3, 4, 5, there hold partial fraction expansions 

1 (5+M)Hfe-5H,.-fc-M-f/„+fc 



(^)n+l 



A:=0 

1/2 
(a;-|-fc)3 

1/6 




(a; + kY 



1/24 



{x+ky 
{ 



x + k 



+3 { i5+t^)Hk -5H„_fc + } X { (5+M)Hf > +5fff/ ^ } 



(5 + M)fffc-5ff„_fc-nH„ + fc} +3{(5 + f)«"i^*+5«"i^ifc-M-H^^_|?t.| +6(5 + Ai)ff 



+ 6{(5 + f.)Hfc-5ff„ 



and the six corresponding harmonic number identities: 
.5 / „ 4- k\t^ 



k=0 



+8{ (5+M)-H"fc-5H„_fc-MH,^+fc} X {(5+^')-f^f ^ -5fff/;,-Atfff_i? J -6Mffl!|!fc 



By means of the standard partial fraction method, we can also derive the following 
algebraic identities and the corresponding harmonic number formulae, even though 
they are not consequences of the theorem and the corollary proved in the present 
paper. 

Example 6. Partial fraction decomposition formula 
x{l - x)"^ 1 / n \^ /■ n + k\^ 



x + T.[k) ( k ) 



-k ^ 1 + 2fcgn+fc + 2fcg„_fc - 4gfc ) 
(2; + fc)2 x + k J 

and the corresponding harmonic number identity associated with Beukers' conjecture 
(cf. P and El); 

= f^(^)'(" + ^)'{l + 2fcff„+fe + 2fcF„_fc-4i?fe|. 

k=l 

Example 7. Partial fraction decomposition formula 

ii+x)i ty) \ k ) 

and the corresponding harmonic number identity: 

n(n + l) f]fc^(^)'("J^)'{i+i/„+,+//„_,-2/7.}. 

fe=i 



8 



CHU Wenchang: Partial-Fraction Decompositions and Harmonic Number Identities 



Example 8. Partial fraction expansion with denominator polynomial of type "2+1" 



n! X (2n)! ^ " C^J (-1)^ 

+ Z](fc)(n + fc){(^ 

and the corresponding harmonic number identity: 



k)' 



+ 



x + k 



^ 1 + n f"+'=) 

fc=l \n+lJ k=0 



n 



2n 
n + k 



{2Hn-k — Hk — Hn+k)- 



Example 9. Partial fraction expansion with denominator polynomial of type "3+1" 



1 



n / 2n 



C+J (-1)'= , ^ 



X (2n)! 

~ (l+n)^ ^ ^n+fc^^ 



+ 



fe=0 



2 / 2n 
n + A; 



+ (.x + fc)2 x + k 

and the corresponding harmonic number identity: 



^ ^ ^ fe=l i„+ij fc=0 



(^2Hk + Hn+k — ^Hn-k) 

n + k J\ + (2Hi'^+Hi%+3Hl%) 



Example 10. Partial fraction decomposition formula with denominator polynomial 
of type "3 + 2" 



n! X {(2n)!}^ 



1 " 

— E 

4- n i-^ 



\n+k) 



{x)l^,{l-x)l 1 + (III) \{x-k) 



fc=0 



X — k 



2n \2 
n + k) \{x + kf 



{x + kf 



+ 



1 



2{x + k) 

and the corresponding harmonic number identity: 



(2) 



■)]} 



n / 2n \ 
^ N ^ Vn+fcj 



E J^n+k\ i^^n+k — Hk-1 — 2iJ„_fc). 



Example 11. For nonnegative integer with < < 4 + 4n, there hold partial 
fraction decomposition formulae: 



(n!)V 



|4 

'"+1 fe=0 



k 

(-fc)«-2 



^(-fc) 



0-1 



e - 4k{Hk - Hn-k) 



{x + kf {x + kf 

2{x + kY L^-^'=(^*-^"-'=^}'-^-"''(^^"+^--)} 

(-fc)^-3 r|e^fc(g,-g„_,)|V2|9-4fc^^(Hf>-gfij| 

6(a; + /c) [ -3{e^fc(i/fc-i/„_fc)}{e^fc^(fff' 



(21a) 

(21b) 
(21c) 
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The corresponding harmonic identities read as 

, 4 



k=0 



{ e-4fe(fffc -_f/„_fc ) } +2{ e-4fe^ (H<''>-_ff^'i J } 
-3 { e-4fc(ff fc - H„ _ ) } { e-4fe" (i/f J } 



0, 



0<9 <2 
6' = 3- 



- 4n 

4?!. 



(22a) 



(22b) 



For 9 = 0,1, 2, the corresponding results to this identity have been conjectured by 
Weideman |10[ Eq21] and confirmed by Driver et al ^ Eq20]. In particular, we 
recover, with the case = 1, the identity found by Driver et al |S1 Eq21]. 



The list can be endless. However, we are not bothered to extend it further. The 
interested reader can do that for enjoyment. 
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Appendix A. Table for JI-Coefficients Computed via ltTT|l and ifT?^ 



r!o(A,/i,-fc) = 1. (Al) 

r!i(A, ^JL, -k) = Aji/fc - Hn-k] + ^i{Hk - Hn+k\- (A2) 

n2{X,^x,~k) = {\{Hk-Hn-k)+^i{Hk-H,,+k)]^ (A3a) 

+ A{i/f (A3b) 

n^{X,^x,-k) = i^X{Hk - H„^k) + It{Hk - Hn+k)} (A4a) 

+ 2{ A(i/f > - Hi%) + > - Hi%) } (A4b) 

+ 3[x{Hk - Hn-k) + i^{Hk - Hn+h) } (A4c) 



10 CHU Wcnchang: Partial-Fraction Decompositions and Harmonic Number Identities 



ni{X,fi,~k) = {x{Hk - Hn-k) + t^{Hk - Hn+k)} (A5a) 

+ 6{x{Hjt^ + Hi%) + ,{Hlt' - } (A5b) 

+ s[\{Hu - Hn-k) + ^i{Hk - Hn+k) } (A5c) 

x{A(i/f -i/fi,)+M(i^f -i^ll)} (A5d) 

+ 6[x{Hk - Hn-k) + ^i{Hk - Hn+k) y (A5e) 

x{A(i/f (A5f) 

+ 3{A(i7f ^ + H<:\) + ) - <p y. (A5g) 

n5i\fi,-k) = {x{Hk - H,,^k) + f^Hk - H„+k)} (A6a) 

+ 24{ A(ijf > - H^%) + ^{Hi'^ H^l) } (A6b) 

+ lQ[x{Hk - Hn-k) + f^{Hk - (A6c) 
x{A«+i/,fJj+Mi/f (A6d) 

+ 20[x{Hk - Hn-k) + fi{Hk - Hn+k)y (A6e) 

X {x{Hi'^ H^\) + f.{Hj^^ } (A6f) 

+ 15{x{Hk - Hn-k) + f^{Hk - Hn+k) } (A6g) 

x{A(i?f+i/,fi,)+/.(Hf (A6h) 

+ 30{x{Hk - Hn-k) + - Hn+k) } (A6i) 

x{A«+H,ri,)+M«-<,)} (A6j) 

+ 20{A(i/f ) + + ^(i7f ^ - Hi%) } (A6k) 

X { A(i7f > - Hi%) + ^,{Hi'^ - Hi%)}. (A61) 



Appendix B. Table for nj-CoEFFiciENTS Computed via H2()a|i 

wo{X,-k) = l. (Bl) 

wiiX,-k) = x{Hk-Hn-k}- (B2) 

W2iX, -k) = X^iHk - + A{i/f + h\%]. (B3) 

W3(A, -k) = A^l^fe - i/„_fe|' + 2x[Hf^ - i/fi^} (B4a) 

+ iX^[Hk Hn-k\ X ' + H^n-k]- (B4b) 

^c74(A, -fc) = X^^Hk - i?„-fc} + 6A{i7f ^ + (B5a) 

+ 8A2{i/fe - H„_,} X {ijf ' - (B5b) 

+ 6A3{i/fe - Hn-k]' X {i/f + i/fi,} (B5c) 
+ 3A^{i/f (B5d) 
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n75(A, -k) = X^^Hk - Hn-kY + 24A{i/f ' - Hi%} (B6a) 

+ 10X^{Hk - X {i/f ) + hI^\^ (B6b) 

+ 2QX^{Hk - }' X {ijf ^ - hI^\} (B6c) 

+ 15X'{h, - i/„_, } X {i/f ' + Hl^\y (B6d) 

+ 30A2 {iJfc - Hr,-k } X {i/f ' + H^;"!^ } (B6e) 

+ 20A^{i7f +i7(^2,}x{i/f -i^i^i,}. (B6f) 

Appendix C. Table for ^-Coefficients Computed via (|20b|) 

uoif^, -k) = 1. (CI) 

wi(/i,-/c) = /i{i/fc-i/„+fc}. (C2) 

u;2(m, -A;) = f^^{Hk - Hn+kV + > - H^^l,}. (C3) 

u;3(m, -k) - /i'^fc - ^^n+fc| + 2/i{i/f ^ - i^fi,} (C4a) 

+ 3M^{i/, - if„+fe} X {i/f > - J. (C4b) 

a;4(M,-fc) = M^{^ffe-if„+fe} + Qf^{Ht^ - H^lk} (C5a) 

+ 8m' { H, - H„+k } X { iff > - } (C5b) 

+ 6^^'{H, - H„+fe}' X {i/f > - (C5c) 

+ 3M'{i/f (C5d) 

W5(/^, -A:) = ^i^{Hk - Hn+k} + 24Ai{i?f ^ - Hj^%} (C6a) 

+ lO/i^ji?, - Hn+k}' X {i/f > - H^^l,} (C6b) 

+ 20//|i/;, - Hn+k}' X {iff > - i//;^ J (C6c) 

+ 15fj,^{Hk - Hn+k} X {i/f > - (C6d) 

+ 30fi^{Hk - Hn+k] X {i/f > - i?^"^,} (C6e) 
+ 20M'{<-<,}x{i/f (C6f) 
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